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A local interaction between photons can be engineered by coupling a nonlinear system to a transmission 
line. The required high impedance transmission line can be conveniently formed from a chain of Josephson 
junctions. The nonlinearity is generated by side-coupling this chain to a Cooper pair box. We propose to probe 
the resulting photon-photon interactions via their effect on the current-voltage characteristic of a voltage-biased 
Josephson junction connected to the transmission line. Considering the Cooper pair box to be in the weakly 
anharmonic regime, we find that the dc current through the probe junction yields features around the voltages 
2eV = nhcOs, where (Os is the plasma frequency of the superconducting circuit. The features at w > 2 are a direct 
signature of the photon-photon interaction in the system. 

PACS numbers: 74.50+r, 74.81.Fa, 85.25.-j 


I. INTRODUCTION 

Creating strong light-matter interaction attracts increasing 
attention due to both fundamental reasons [1-^] and its poten¬ 
tial application in quantum communication science [1, 5, 6]. 
A prototypical system for studying that interaction consists of 
a quantum system inside a photonic cavity [7]. However, re¬ 
cent rapid experimental advances in several areas [8-13] have 
focused attention on one-dimensional systems in which the 
quantum system is embedded in a waveguide or transmission 
line. In the absence of coupling, photons propagate freely 
down the line. A coupling between the quantum system and 
the line generates an effective photon-photon interaction that 
causes correlations among the photons. This has lead to, for 
instance, the prediction of Kondo physics [14], anti-bunching 
resulting from a photon-blockade effect [15, 16], inelastic 
photon scattering [15, 17, 18], giant Kerr nonlinearities [19], 
and entanglement among photons of different frequencies in 
the line [20]. 

The strength of the coupling between the local quantum 
system and the transmission line has been studied theoreti¬ 
cally in detail in the ohmic spin-boson model, which consists 
of a single two-level system (the spin) bilinearly coupled to 
the photons in the line (the bosons). It was shown [21, 22] 
that the coupling parameter is set by the ratio of the line 
impedance, Z, to the quantum of resistance, RQ = hl{2e)^ « 
6.45k^l. The impedance of typical transmission lines is of 
order the vacuum impedance, Zyac ~ 377 ^1, thereby allowing 
only weak coupling. 

Superconducting circuits are a promising platform for ex¬ 
ploring strong coupling phenomena, and, indeed, the first ex¬ 
periments observing such phenomena have appeared [23, 24]. 
One benefit of using superconducting circuits is that a chain 
of Josephson-coupled superconducting islands acts as a trans¬ 
mission line with a large tunable impedance Z<Rq, which is 
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only limited by the superconductor/insulator transition [25]. 
Recent experiments have studied the microwave properties of 
such Josephson junction chains [26-29]. Moreover, supercon¬ 
ducting circuits allow the realization of a variety of quantum 
systems that behave like artificial atoms [1, 30]. 

In our work, we take the quantum system coupled to the 
transmission line to be a Cooper pair box [31, 32]. Then, we 
propose to detect the photon-photon interaction generated by 
that system by measuring the dc current-voltage characteristic 
of an additional Josephson junction connected to the transmis¬ 
sion line. According to dynamical Coulomb blockade theory 
[also called F(£’)-theory] [33], Cooper pairs can tunnel inco¬ 
herently through that probe junction provided that they can 
release their energy 2^y into the environment, which in our 
case consists of the transmission line with the side-coupled 
circuit. Therefore, the dc current refiects both the elastic and 
inelastic scattering properties of photons. 

Let us consider the current-voltage characteristic in more 
detail: In the harmonic regime, the effective impedance of 
the environment is almost fiat, except at frequencies near the 
plasma frequency COs of the superconducting side circuit (the 
Cooper pair box). This results in a feature at 2^y = hcOs in the 
current-voltage characteristic. Anharmonic corrections cause 
additional features near the voltages 2^y = nhcOs {n>2 in¬ 
teger) through inelastic photon scattering. These additional 
features are most pronounced in the strong coupling regime 
when Z approaches Rg (Z < Rg). 

The paper is organized as follows. In Sec. II, we intro¬ 
duce the Hamiltonian that describes the circuit studied. The 
current-voltage characteristic of the probe junction in the har¬ 
monic regime is calculated in Sec. III. In Sec. IV, we include 
a weak anharmonic correction and study the effect of the re¬ 
sulting photon-photon interaction on the current. Finally, we 
conclude in Sec.V. 
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FIG. 1. The system consists of a transmission line that is capaci- 
tively coupled (capacitor Q) to a Josephson junction, shown inside 
the dashed box. The transmission line is realized using a chain of 
Josephson junctions with Josephson energy Ej much larger than the 
charging energy Ec- The system is probed at node m using another 
Josephson junction (outside the dashed box) whose current-voltage 
characteristic is sensitive to the properties of the system. 


II. THE CIRCUIT STUDIED 


We are interested in the interactions of photons propagating 
in a non-linear electromagnetic environment. In particular, we 
study a transmission line, consisting of a chain of Josephson 
junctions, to which an additional Josephson junction acting as 
the non-linear element is side-coupled at node ^ = 0 as shown 
in Fig. 1 (dashed box). We assume weak coupling, namely the 
coupling capacitance, Q, is much smaller than the character¬ 
istic capacitances of the chain and the non-linear element. 

The Hamiltonian of the system is, thus, assembled from 
three parts, 

H = (1) 


where Hj is the Hamiltonian of the transmission line, Hj is 
the Hamiltonian of the side-coupled Josephson junction, and 
He is the coupling Hamiltonian. 

A transmission line with large impedance can be realized 
using a chain of Josephson junctions in the limit where the 
Josephson energy is much larger than the charging energy 
Ec [26, 27] . The chain is described by the charge and phase 
operators at each node n, denoted Qn and respectively. 
They are conjugate variables satisfying the commutation rela¬ 
tion [Qn-,^m\ = —2ie5nm- As Ej ^ Ec, phase fluctuations are 
small and we may approximate the Josephson coupling by a 
quadratic term. We further consider the case where the capac¬ 
itance to the ground Co is much larger than the mutual capac¬ 
itance C. Then, for frequencies much smaller than the plasma 
frequency of Josephson junctions in the chain, the Hamilto¬ 
nian takes the simple form [25] 


Ht=E 


2Co 


+ 


1 {^n 

( 2^)2 2 Z 


( 2 ) 


where the inductance is L = 1/ {Ae^Ej). Note that we use units 
where fi=\. At frequencies m <C (5)o, where COq = I/^/LCq, 
the transmission line has a linear spectrum. 

The side-coupled Josephson junction with Josephson en¬ 
ergy Ej is described by the Hamiltonian 


(Qs+CgVg)^ 

ELj = 


2Ce 


— Ej cos 05 


( 3 ) 


where and 05 are the conjugate charge and phase operators 
at node 5 (see Fig. 1). Furthermore, Cg and Vg are the gate ca¬ 
pacitance and gate voltage, respectively, and Cj: = + Cg is 

the total capacitance of the side-coupled Josephson junction. 

Finally, we turn to the coupling Hamiltonian He. When 
the coupling capacitance is small, Q <C Co,C 5 :, the coupling 
Hamiltonian reads 


He = 


Cc 

CqCy. 


Qo(Qd 


(4) 


where we used the fact that for C <C Q the coupling is local, 
i.e., the side-coupled Josephson junction couples only to the 
charge Qo Sit n = 0. The Hamiltonian H fully describes our 
non-linear system. 

As a next step, we introduce the probe circuit used to 
characterize the photon-photon interactions generated by the 
non-linear system. The probe circuit consists of yet another 
Josephson junction, with Josephson energy Ej and in series 
with a voltage source as shown in Fig. 1, coupled to the trans¬ 
mission line at node m [34, 35]. The current-voltage charac¬ 
teristic of the probe Josephson junction is influenced by the 
correlations of the phase 0 m (0 ^^de m, correlations that de¬ 
pend on the fluctuations in the non-linear environment. The I- 
V characteristic may, thus, be used to characterize the photon- 
photon interactions in the non-linear system. 

In particular, using /’(E')-theory, it can be shown that at zero 
temperature the current flowing through the probe Josephson 
junction takes the form [33] 


l{y)=ne{E]fp{2eV), (5) 


for voltages eV < 2A, where A is the superconducting gap, and 


P{E) = E [dt (6) 

—oo 


is the probability of the probe Josephson junction to emit en¬ 
ergy E to its environment, described by the Hamiltonian F/env 
Though -theory is usually presented in the context of a 
linear environment (F/env is assumed to be quadratic) [33, 36- 
38], Eqs. (5)-(6) hold more generally for a non-linear environ¬ 
ment [34]. If the capacitance of the probe Josephson junction 
is sufficiently small, ^ Co, the Hamiltonian of the environ¬ 
ment in Eq. ( 6 ) may be replaced by the Hamiltonian H of the 
non-linear system we want to characterize. Our task is then to 
compute the phase correlator 


III. THE LINEAR REGIME 

As a first step, we will consider the system in the linear 
regime, where photons do not interact. That is, we assume 
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Ej ^ j {2Cz) and approximate the junction Hamiltonian Hj 
in Eq. (3) by 


„(o) _ {Qd+CgVgf 
-2ft 


2 




(7) 


In this section, we study the behavior of this simplified system 
described by = Hr + He to set the basis for in¬ 

vestigating interaction effects, the main focus of our work, in 
the following section. In this regime, the gate voltage Vg can 
be gauged out of the Hamiltonian, and the side-coupled cir¬ 
cuit behaves as an harmonic oscillator with plasma frequency 
(Os = 2e^E]jCz- We will assume that (Os ^ (Oq. 


A. Phase-phase correlator 


As the system is non-interacting, the phase-phase correla¬ 
tion function in Eq. (6) can be simplified by exploiting Wick’s 
theorem [33]: 

( 8 ) 

where 

•7(0 = (9) 

To evaluate the correlator, we use the retarded Green func- 

tion = i@{t){[<l)n{t),<l>m{ 0 )])H( 0 )^ where 0(r) is 

the Heaviside step function. The relation between the two is 
most easily written in frequency space. At zero temperature, 
it reads 


(<^>„(O0„(O))«(o)=2|ge- 




. ( 10 ) 


The local Green function (o) needed to com¬ 

pute J{t) is obtained by deriving its equation of motion and 
using scattering theory (see Appendix A). Photons propagate 
freely in the transmission line and are scattered by the side- 
coupled harmonic oscillator at node n = 0, yielding the refiec- 
tion coefficient 


r{(0) 


l-2i^ 

(0 



CoC^ (iP- - (Up 



In terms of this refiection coefficient, the Green function is 


^ (0m 5 0m? to) 


,71 Zo 
( ORq 


\^r{(o)e ^ 


( 12 ) 


where Zq = ^JlICq is the impedance of the chain. Under the 
conditions specified above, C^CzjC^ '> 1 and (Os ^ (Oq, the 
refiection coefficient has a narrow resonance at (0 = (Os with 
width 


r 


1 eg COs 

4 CqCz (Oq 


(Os. 


(13) 


Close to the resonance, we can approximate Eq. (11) as 
r{(0) = -l/[l-i{(0-(0s)/r]. 
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FIG. 2. (Color online) The linear regime: Current-voltage char¬ 
acteristic of the probe Josephson junction when placed at m = 0. 
The parameters are Y/(Os = 0.02, Eeui-oii/(Os = 20, and different 
Zq {ZqIRq = 0.01, 0.1, 0.2). The side-coupled Josephson junction 
causes a resonance at 2eV = ft)^. In the limit ZqIRq^ 0, the current 
vanishes at the resonance. 


Substituting Eq. (10) into Eq. (9) and using the Green func¬ 
tion (12), one obtains 


■7(0 = |-/^9?[Z(a))](e-'«-l), (14) 

^ 0 


as expected from /^(E")-theory [33], with the impedance 


z{(o) = 


2 


\-Yr{(o)e ^ 


(15) 


The prefactor 1/2 corresponds to the fact that the probe junc¬ 
tion ‘sees’ an environment consisting of two half-infinite 
transmission lines. Ear from the resonance at (Os, the 
impedance is unaffected by the side-coupled Josephson junc¬ 
tion as r{(o) 0. In contrast, at the resonance, photons are 

strongly scattered. In particular when the probe and the scat- 
terer are coupled to the same node (m = 0), r{(Os) = —I so 
that transport is completely blocked due to destructive inter¬ 
ference. Changing the distance between the probe and the 
scatterer modulates the phase difference between incoming 
and reflected photons and, thus, creates an interference pat¬ 
tern. 


B. Current-voltage characteristic 

To compute the current-voltage characteristic, we need to 
determine P{E). This can be done numerically using the inte¬ 
gral equation [33] 


n, 

EP{E) = 2 jdE 


,^[Z{E-E') 


Rq 


P{E'). 


(16) 


The result can be obtained by starting with an arbitrary 
value P(0) and then using the normalization condition 
/o^cu'-off= I, where ^cut-off > (Os. 

The current-voltage characteristic is plotted in Eig. 2 for 
several values of the impedance of the transmission line. The 
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FIG. 3. (Color online) The resonance in the current-voltage charac¬ 
teristic for different values of the distance between the side-coupled 
Josephson junction and the probe Josephson junction, a = 2m(Osl (Uo- 
Results in the single-photon, linear regime are plotted [Eq. (18)] with 
T/(Os = 0.02. Note the effect of interference on the shape of the res¬ 
onance. 


characteristic current is given by 


ne{E]f Zq 
(Os Rq' 


(17) 


The background current decreases with increasing voltage. In 
addition, there is a clear resonance feature at 2eV = (Os. 

This result can be understood as follows. The starting point 
is to recognize that when a bias voltage V is applied, Cooper 
pairs can flow through the probe junction provided that they 
can release their energy by emitting one or several photons 
into the environment. 

First, let us concentrate on the regime Zq/Rq<^\. In that 
case, multi-photon processes are suppressed, and we can ex¬ 
pand 1 +/(^). Thus, the current is proportional to the 

Fourier transform of J{t) at frequency 2eV. It is straightfor¬ 
ward to show that for a constant impedance this yields a cur¬ 
rent that decays with increasing voltage as /(V) oc I/y. On 
top of this, the resonance in the impedance at (Os due to the 
side-coupled Josephson junction leads to a resonance in the 
current-voltage characteristic at 2^y = (Os. Namely, the cor¬ 
rection to the current 5/i (SV) at voltages V = (0st(2e) + SV 
takes the form 


5/i(5y) 


-r^ 


(2^5y)2+r2 


2^5y . . 

cos a-z^sina , (18) 


where a = 2m(0slc0o. This leads to a complete extinction of 
the current at 5y = 0 (at the one photon level) when the probe 
is coupled to the same node as the side-coupled Josephson 
junction (m = 0). The shape of the resonance for different 
a is shown in Fig. 3; note the sensitivity to the placement of 
the probe produced by interference effects. The width of the 
resonance is given by W\ = r/e where F is given in Eq. (13). 


Let us now turn to multi-photon processes corresponding to 
higher order terms in J{t). These processes modify the reso¬ 
nance at 2^y = (Os. In particular, while the scattering from the 
side-coupled Josephson junction may completely block the 
single-photon process at that voltage, this is not the case for 
the multi-photon processes: at most one photon can be on res¬ 
onance, whereas the other photons will be off resonance and 
therefore propagate freely. Thus, the multi-photon processes 
lead to a finite current at the resonance. As an ^-photon pro¬ 
cess yields a current contribution proportional to {ZoIRqY, 
the resonant structure weakens with increasing ZqIRq due to 
the increasing importance of multi-photon processes. 

In addition, one might expect that multi-photon processes 
lead to higher order resonances at voltages 2^y = n(Os {n>2). 
We find, however, that this is not the case. While 2^y = n(Os 
is indeed a resonance condition for an ^-photon process, the 
non-resonant background from the entire frequency range is 
large enough to completely overwhelm that contribution. 

Thus, in the linear regime where photons do not interact, the 
side-coupled Josephson junction leads to a single resonance 
in the current voltage characteristic at 2^y = (Os. As we will 
show below, additional features at 2^y = n(Os with ^ > 2 are 
a signature of photon-photon interactions. 


IV. THE NON-LINEAR REGIME 


To investigate photon-photon interactions, we now take into 
account the non-linearity of the side-coupled Josephson junc¬ 
tion. In particular, we concentrate on the case of weak non¬ 
linearity in the regime j (2 Cy). To do so, we expand 

Eq. (3) up to fourth order in 0^, 

(19) 

where 

= (20) 
In the following, we treat y as a perturbation. 


A. Phase-phase correlator 

As the Hamiltonian describes an interact¬ 

ing system, we can no longer write a closed form expression 
for the phase-phase correlator (6) in terms of (0n(O^m(O))- 
Instead we expand (6) in powers of (j)m as follows. 


= 1 + iiut) - uo)]uo))h^,+\{{€{ t) - 
-I +aut) - ^(<?>®). 


( 21 ) 
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FIG. 4. Dyson equation for the two-point Green function. The non¬ 
linearity results in a self energy, E = jl. 


Here the two-point phase-phase correlator represents single 
photon processes, whereas the four-point phase-phase corre¬ 
lators represent two photon processes. As before, we will use 
Green functions to evaluate the correlators. In addition to the 
two-point Green function 0m; we now also need 

the four-point Green functions 0^; to), 0^; to), 

and G/?(0;^,0^;t)L)). In order to facilitate doing perturbation 
theory in the interaction V, we switch to imaginary-time- 
ordered or Matsubara Green functions, ^. 

Let us first evaluate the two-point Green function 
^[0 „(t) 0^(O)], corresponding to single photon processes. As 
shown in Appendix B, using the Dyson equation, we can 
sum up the perturbation series to all orders in the interac¬ 
tions. The representation in terms of Feynman diagrams is 
shown in Fig. 4. For the corresponding retarded Green func¬ 
tion G/?(0;^, 0;^; to), we finally obtain 

0m; ~ ^ 0m; (22) 

^ ^ 0^; ^ (0^5 0m; 


where G^^^ is the Green function in the absence of interac¬ 
tions. 

Using the Green functions G^^^ and (0|)^(o) derived in Ap¬ 
pendix A, we find that the local Green function preserves its 
form though with a shifted resonance frequency m'. Namely, 


Gr (0^,0^; m) 


,n Zo 
'coRq 


n ^ 2/^'” 


(23) 


where 




1-2/^ 


CoCi: - (0, 

Cl co^ 


f2" 


1 -1 


(24) 


with to' (^^[l — (0^/(8£’))]. In the same way, we can show 
that this is true for all two-point Green functions. Note that 
5o)s = (Os/i^Ej) ^ COs coincides with the shift of the excita¬ 
tion energy between the ground and first excited states of the 
Hamiltonian (19). 

Next we turn to the four-point Green functions, correspond¬ 
ing to two-photon processes. Using perturbation theory, we 
may express them in terms of the two-point Green functions. 
As we saw above, it is essential to sum up the perturbation 
series to all orders in V to obtain these two-point Green func¬ 
tions. By contrast, we will keep only the lowest order term in 
V accounting for interactions between the two photons. Then, 
the four-point Green function G/? (0^, 0^; m) has two contri¬ 
butions: the first one describes the independent propagation of 
the two photons, whereas the second one describes the interac¬ 
tion effects. More precisely, the imaginary-time-ordered four- 
point Green function may be written as ^ [0 ^(t)0^(O)] = 

with 


^EjJ dx' [t5iT')tm{0)] • (25) 

0 


The corresponding Feynman diagram is shown in Fig. 5(a). The expression for the local retarded Green function at zero temper¬ 
ature reads (see Appendix C) 

oo 

^ / dcoi 3 [G/?(0^,05;C)l)i)]G/?(0^,05;C)l)±(Oi) 

^0 

The leading order term for the other four-point Green functions Gr [(l)^, Co) and Gr ((l)n, <1)^1 0)) is linear in Ej. In particular, 
we find the local Green functions 





E^ f f 

GR{^li,^m\(^) — dCOi / JC02 3 [G/?(0m,05;^^l)]3 [G/?(0m,05;^^)] ^ G/?(0m, 0^; tO + COi + 5'2t02) (27) 

^ { i ^1.^2=± 


and GR{^rn^^m\(^) = GR((l)^,(l)m',Co). The Feynman diagrams With the above results we can now write the phase-phase 
for the corresponding time-ordered Green functions are shown correlator needed to compute P{E) in the following form 

in Fig. 5(b). 


^ /it) ^ 


(28) 
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with 


J\t) = - f d(o2[GRi(l)„„(l)„,-,co)] (e-'® -1) (29) 

71 J 
0 
oo 

■“‘(r) ~ 11^(0 { is [S&M, «)] (30) 

0 


B. Current-voltage characteristic 


Using Eq. (28) to compute P{E), we obtain the current 


I{V) ~ e {E^f {\ J [i2eVl + j'{l)] (31) 

'8Gf\(^l(^l-2eV)] - is [G«(<^)3,,^)„;2ey)] }. 




The first line describes the resonant structure discussed in 
Sec. III. Here the only effect of the non-linearity is to shift 
the resonance from (Os to to'. The second line describes inter¬ 
action effects between two photons. The current-voltage char¬ 
acteristic including these effects is shown in Fig. 6: it displays 
additional structure at 2^y = 2(0'. 

The new peak at 2^y = 2(o' comes from the contribution 
^ 5G^^(0^,0^;2^y). This contribution describes a process 
in which a Cooper pair tunnels through the probe Josephson 
junction emitting two photons. When both photons are on res¬ 
onance with the side-coupled Josephson junction, they inter¬ 
act strongly. This happens when each photon takes away half 
of the energy of the Cooper pair, (O = ^y (o'. The resulting 
correction to the current is obtained using Eq. (26). As shown 
in Appendix C, for voltages V = (o'Je-\- 5V, it takes the form 


5l2{5V) = 


X 


n Zo co', 

^ 32 RqE} [{e5Vf + r^Y 

Ye5Vcos{2a')-]^[{e8vf 


(32) 


r^] sin(2a') 



FIG. 5. The Feynman diagrams for the interaction correc¬ 
tion to the four-point Green functions, (a) [0 ^(t)0^(O)] . 

(b) ^ [(l>^{T)(l>m{0)] and ^ [(l>^{T)(l>m{0)]. 


where a' = 2m(o'JcOo. 

The characteristic amplitude A 2 of the change in current 
is, thus, much smaller than Iq or the single-photon resonant 
structure 5/i, 


A 2 


71 Zq CO'^ 


/o</o. 


(33) 


Here, the suppression factor Zq /Rq is due to the fact that it is 
a two-photon process, whereas the suppression factor (o'jEj 
is due to the fact that it is an interaction effect. Notice that 
the widths of the resonances at 2^y = (o' and 2^y = 2(o' are 
the same. The dependence of the shape of the second reso¬ 
nance on the distance oc a' between the side-coupled Joseph¬ 
son junction and the probe Josephson junction is shown in 
Fig. 7. 

We finally consider the current contribution stemming from 
;(()). While it is in principle of the same order as 
the current contribution from 5G^^(0^, 0^; (o), i.e., it is pro¬ 
portional to (Zo//?g)((o'/£’))/o, in this case it is impossible to 
fulfill the resonance condition simultaneously for all the pho¬ 
tons involved. Therefore, this contribution acquires an addi¬ 
tional suppression factor E/ol)', and we can neglect it. 

The main interaction effect is, thus, the appearance of a res¬ 
onance at 2^y = 2(0' due to two-photon processes. Higher 
order processes are expected to lead to additional features 
at 2^y = n(o'^{n > 3). However, their amplitude decreases 
rapidly with increasing n and may be estimated as An ^ 
[(Zo/«e)(coV£))]”-‘/o« A2. 


V. CONCLUSION 

We have shown that the dc current-voltage characteristic 
of a Josephson junction provides a sensitive probe to study 
photon-photon interactions in a non-linear environment. In 
particular, we investigated the case of a transmission line side- 
coupled to another Josephson junction whose non-linearity 
leads to local photon-photon interactions. Scattering of in¬ 
dividual photons by the side-coupled Josephson junction re- 



FIG. 6. The non-linear regime: Current-voltage characteristic of 
the probe Josephson junction when placed at m = 0. The parame¬ 
ters are cd'/E) = 0.9, T/ft)' = 0.02, = 20, and ZqIRq = 

0.2. Photon-photon interactions lead to a second resonant feature 
at 2^y = 2 ft)'. A zoom on that feature with amplitude SI 2 /I 0 
(Zq/Rq){(o'JE^) is shown in the inset. 
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FIG. 7. The second resonance in the current-voltage characteristic 
for different values of the distance between the side-coupled Joseph- 
son junction and the probe Josephson junction, a' = 2m(o'J(OQ. Re¬ 
sults are plotted near leV = 2ft)' in the two-photon, non-linear regime 
[Eq. (32)] with F/m' = 0.02. 


see how these features are modified in the strongly non-linear 
regime. 


ACKNOWLEDGMENTS 

The work of Li-Jing Jin is supported by the China Scholar¬ 
ship Council. MH and JSM acknowledge support by ANR 
through grants ANR-11-JS04-003-01 and ANR-12-BS04- 
0016-03, and an EU-FP7 Marie Curie IRG. HUB acknowl¬ 
edges support from U.S. NSF Grant No. PHY-14-04125. FH 
is supported by Institut universitaire de France and by the Eu¬ 
ropean Research Council (grant no. 306731). We thank the 
Fondation Nanosiences de Grenoble for facilitating the ex¬ 
change between Grenoble and Duke. 


suits in a resonance feature in the current-voltage character¬ 
istic of the probe Josephson junction at 2^y = ft)', where ft)' 
is the plasma frequency of the side-coupled Josephson junc¬ 
tion. By contrast, the interactions due to the non-linearity 
yield an additional resonance feature at 2^y = 2 ft)' due to two- 
photon processes. Such a feature is thus a clear indication of 
photon-photon interactions. While we concentrated here on 
the regime of a weak non-linearity, it will be interesting to 


Appendix A: Two point retarded Green functions 

In this appendix, we derive the two-point retarded Green 
functions of the linear system using equations of motion 
and scattering theory. 

We start with the coupled equations of motion for 
{(j>n,<l>m', (O) and (<p5,<pm; CO), 


((/»«,<?>,«; 0)) - co^ (1 + ((/»„,(/»„; w) - Gf^ ((/»„+i, 0™; (o) - G^^ o)) 

- -p^ (<^>5; ®) ^nO + ^ ^nm = 0, 

(-0 (-0 


- (O^) Gf' ((t>S,^m-, ffl) - CO) - G^ ((^>1, “ G^ «) 

Ce l 


= 0 . 


(Al) 


(A2) 


Combining Eqs. (Al) and (A2) then yields the following equation for G^^^ (<pn, (/>„; (O), 


C0^G‘'^\(j>n,(l)m;C0)-C0^[ 1 + 


ci 


(0^ 


CqCz 


'^nO 


2Gjj”^ ((l)n, (O)- (</»«+!, <l>m; (0)-G^^^ (O) 


Co 


5™. (A3) 


If there is no side-coupling, Q = 0, Eq. (A3) describes pho¬ 
tons propagating freely along the infinite transmission line. At 
frequencies ft) ^ ft)o, the dispersion is linear, ft) = cOok with 
wavevector ^ ^ 1, and the solution is 


r(ft)), as well as the amplitudes A and B have to be determined 
using the boundary conditions at ^ = 0 and n = m. One finds 
t{co) = 1 -\-r{co) with r{co) given by Eq. (11) in the main text. 
Furthermore, 


^(0) 

chain 


(^ni 0m? to) 


i— ^ik\n—m\ 

CO Rq 


(A4) 


The side-coupling leads to scattering of photons at the node 
n = 0. Then, for m > 0, the solution may be written in the 
form 




n>m, 

< B + r{(0)e^^^] 0<n<m, (A5) 

Bt{co)e-^^^ n<0, 


where the refiection and transmission coefficients, r{co) and 



(A6) 

(ORq 

A=5[e-2'^'” + r((0) . 

(A7) 


The result is obtained by substituting Eqs. (11), (A6), and 
(A7) into Eq. (A5). Generalizing to arbitrary m, we find 


i<Pn, <!>„,■, (0) = i-^ + r(m)e'^(l"l+l'”l) 


(OR, 


(A8) 
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The local Green function needed to evaluate P{E), thus, reads 


(OKq L 


(A9) 


While this is the only Green function needed in the linear 
case, more Green functions are required in the non-linear case. 


Using Eq. (A2), we obtain 

of (AlO) 

UJ axQ^qUJ 

Similarily, the Green functions 

obey coupled equations of motion. One 
may show that (0^, 0^; (o) = whereas 

the equation for G® (0^, 0^; m) is 


{(0^ - (oj) G® (<^>5,05;ft)) - (00,0^;«) - i<t>i ><(> 5 ;®) - {<t>-u<t>5'^0)) 

Ce L 


.M 


(All) 


Using Eq. (AlO), one obtains 


(05,05;a)) = -4i 


0)RqZo{Cc(0) 


rr(co). (A12) 


Using the explicit expresssion for r{co), Eq. (A12) may be 
rewritten as 


(05,05;®) = “1 


1 


RqC^ fi )2 _o )2 + i^^^o)- 2/6)0) 

(A13) 

Finally, using the fact that C^KCqCz) ^ 1, we approximate 


(05:05;^^) - 


271 


RqCz CO^ — {(Os — iE) 


2 ’ 


(A14) 


where 


r= - 


1 C? m 2 


4 CqCe tOo 
This result also allows us to evaluate 
1 /. 

(05)«(O) = [g^^^ (05,05;®) 


( 0 , 

= w 


Appendix B: The Dyson equation 

In the following, we present the derivation of Eq. (22). 
The time-ordered two-point Green function is defined as 


= {T,U^)^m{0))H,,, (Bl) 


where T^; is the time-ordering operator. Eq. (Bl) can be rewrit¬ 
ten as 


^t[0«(t)0,„(O)] 


(7^0^ (t)0^ ( 0 ) 5 ( 00 ) )^(o) 
(‘^H)//(o) 


(B2) 


where 5(oo) = exp [— /q°° Jt' V (t')] . 

Expanding Eq. (B2) up to first order in the perturbation V = 
—Ej(j)^/24 and using Wick’s theorem yields 


[0«(T)0m(O)] [0„(T)0™(O)] + J dx' ) [0„(t)05(t')] y (0|)H(O)^f ^ [05(T')0m(O)] • (B3) 

0 


After Fourier transformation and analytical continuation, one 
obtains the corresponding retarded Green function, 

GR{^n,^m\(0) — G^R\(^n,^m\(0) (B4) 

-f G® {(\)n,(\)5\(0) ^ {(i>8,(\)m\(0). 

While far from the resonance dX CO = (Os the second term is 
much smaller than the first one, this is no longer true close 
to the resonance. Thus, this first order expansion is not suffi¬ 


cient is to describe the modifications to the resonance due to 
the perturbation. It is possible to go beyond the first order ex¬ 
pansion by realizing that £’)(0|)^(o) /2 is a local self-energy, 
0^). Thus, one obtains the Dyson equation 

= GR\(l)n, (l)m', (O) (B5) 

+Gj; (0n, 05 ; CD) y ( 0 I) ^(o) GR{^s,^m',0)) , 

A similar equation can be written for the Green function 
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GR{(p5,<l>m',(0)- Namely, 

GR{^S^4>m',0}) = {([)§,^m',(0) 

^JO) , 


Combining Eqs. (B5) and (B6), we obtain the result 


(B6) 


+Gs ' {<l>d,<l>S'->0)) ^{(l)§)ff( 0 )GR {(l>s, (pm', 0}) ■ 


/- I'A A A i^n,h',(>^)^{(pl)H(fi)^''R ('p5,(pm',(»y) 


E^r 


z(^) I 


(B7) 


Then, using the Green functions of the linear problem 
derived in appendix A, we find that the full Green func¬ 
tion GR{^n^^m\(^) has the same form as G]^\(j)n^(l)m\(o) 
though with a shift of the resonance frequency, (Os ^ co'^ ^ 
(l — (0|)^(o)/4). Similarily, we find that this frequency 
shift appears in all two-point Green functions. 


Appendix C: The contribution of photon-photon interaction 

The four-point retarded Green function 
needed to compute the interaction contribution to the current- 
voltage characteristic, is obtained from Eq. (25) by taking the 
Fourier transform and then performing the analytical contin¬ 
uation from Matsubara to real frequencies, icOy (0 /0+, 
and using standard methods of contour integration. One 
obtains Eq. (26) which takes the form 5G^^ (0^, 0^; m) 


{Ej/7r^)f^{(o), where 


/(^) = 1 ^ dCOiZ 

(Cl) 

The integral is dominated by frequencies where both Green 
functions are close to resonance, (Oi ^ (O ± (Oi ^ co'^. This 
requires (0 ^ 2(o'. We, thus, approximate (0 = l(o[ -f 5m and 
mi = m' + 5mi. The Green functions close to resonance take 
the form 




1 


1 


(o'^RqCcCo'.I-i^ 




We then rewrite 


/(«)- 


It: 


RqCc{(o's)^ 


i^\m\ 


7 sin^!^ + ^cos^!^ 

J dScOi - r «o 


l + (^)‘ 


1-i 


dw-dcoi 

r 


It is straightforward to evaluate the convolution integrals to obtain 


( TtCc 5mcosa'-f2rsina' 


3 [/(m)] 




nCr 


\ IRqCqCy, 

Finally, to compute the current-voltage characteristic, we need 

2E\ 


2T cos a' — 5(0 sin a' 
(5m)2-F4r2 


"I / h ^ 

{<pUm',M + 50))] c. ^9? [f{ 2 ( 0 i + 5o))] 3 [f{2(oi + 5o))] 

=" "¥5 (§) [(5«)^+4P|^ {25<»rcos(2a') - i -4r^)si„(2a')}. 


(C3) 


(C4) 

(C5) 
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